Interactions of Spin Waves with a Magnetic Vortex 
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We have investigated azimuthal spin-wave modes in magnetic vortex structures using time- 
resolved Kerr microscopy. Spatially resolved phase and amplitude spectra of ferromagnetic disks 
with diameters from 5 /im down to 500 nm reveal that the lowest order azimuthal spin wave mode 
splits into a doublet as the disk size decreases. We demonstrate that the splitting is due to the 
coupling between spin waves and the gyrotropic motion of the vortex core. 
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The magnetic ground state of a soft ferromagnetic par- 
ticle is determined by a balance between magnetostatic 
and exchange energies. In the case of disks with negligible 
magnetocrystalline anisotropy, the stable configuration 
for radii of several hundred nanometers and thicknesses 
of the order of 20 - 50 nm is a vortex with a core diameter 
on the order of the exchange length0,|3|. The magnetiza- 
tion outside of the core circulates around the central axis, 
reducing the total magnetostatic energy. This comes at 
the cost of a large exchange energy near the vortex core, 
inside of which the magnetization rotates out of the plane 
of the disk. Because the vortex state represents the sim- 
plest type of domain configuration that can be created in 
a uniformly magnetized particle, its excitation spectra is 
of fundamental interest. Recent experimental work has 
identified two classes of excitations in magnetic vortices. 
First, low-order magnetostatic spin waves can be ob- 
served by either time-resolved Kerr microscopy [3 JEQ 
or Brillouin scattering0,l3- A second type of excitation 
is associated with the translational degree of freedom of 
the vortex core itself 0, 0, E3 

This has a much lower 
characteristic frequency than the other spin- wave modes, 
and it is for this reason that the two types of excitations 
are often treated as distinct. For example, the radial 
spin wave modes of a particle can be calculated accu- 
rately even if the position of the vortex core is regarded 
as axed[l[|. 

In this Letter we report on a study of the dynamics of 
the lowest order azimuthal spin wave modes in ferromag- 
netic disks with diameters between 500 nm and 5 /jm. In 
larger diameter disks, these are degenerate as expected 
for a system with cylindrical symmetry. We demonstrate 
that this degeneracy is lifted in smaller disks (with di- 
ameters less than 2 /zm and thickness to diameter ratio 
> 0.005) due to the motion of the vortex core. The rela- 
tive phases of the two modes are determined by the polar- 
ity of the vortex core, and the magnitude of the splitting 
is of the same order as the vortex gyrotropic frequency. 
These results demonstrate a significant coupling between 
vortex dynamics and magnetostatic spin waves. 

Two sets of permalloy disks were fabricated on Si 
substrates using electron-beam lithography and electron- 



beam evaporation of permalloy (Nio.8iFeo.19). The first 
set of disks were used primarily for investigations of the 
vortex gyrotropic mode and were 50 nm thick with di- 
ameters D from 500 nm to 2 /im. The disks used for 
the detailed investigation of azimuthal spin-wave modes 
covered the range of diameters from 700 nm to 5 /im 
and were 20 nm thick. Each disk exhibited a single- 
vortex ground state as determined by magnetic force mi- 
croscopy. 

Measurements were conducted using time-resolved 
Kerr microscopy 

(TRKM)[HE3- 

Each substrate was 
thinned to a thickness of approximately 25 /im and 
mounted over the 30 /im wide center conductor of a copla- 
nar waveguide, as shown in Fig^a). The waveguide was 
placed on the scanning stage of an oil-immersion micro- 
scope. The TRKM technique probes the response of the 
magnetization to a fast magnetic field pulse h p , which is 
along y in FigOJa). The magnitude of h p is 5 - 10 Oe, and 
the temporal width of the pulse is 120 - 150 psec. The 
polar Kerr rotation 9(x,y;t) of an optical probe pulse is 
measured as a function of the time delay t between h p 
and the subsequent arrival of the probe pulse. A lock- 
in technique is used to measure the change A9(x,y,t) 1 
due to h p , averaged over many pulses. The z-component 
m z (x,y;t) of the dynamic magnetization is proportional 
to Ad(x,y;t). 

In the remanent state, the vortex core is located at 
the center of the disk, as shown in the magnetic force 
micrograph of Fig. ^b). When the pulse h p is ap- 
plied, the core is displaced and then slowly gyrates about 
the center of the particle as the system relaxes towards 
equilibrium 0, 0- This mode, which has also been ob- 
served in a reTonance experiment 0, appears as the 
long-lived low frequency (< 1 GHz) oscillation in Fig. 
Wic), which shows the polar Kerr signal measured at a 
position 250 nm away from the center of two disks with 
diameters of 1.5 and 1 /im. The evolution of the vortex 
mode eigenfrequency as a function of the disk diameter 
is shown in Fig. ^d). As discussed in Ref. fl3l the vortex 
core oscillates in an effective potential that is due to the 
magnetostatic energy of the displaced vortex. The solid 
curve shown in Fig. ^d) follows from a model in which 
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FIG. 1: (a) Schematic of the experiment showing the sample 
on a coplanar waveguide. The current pulse generates a mag- 
netic field pulse along y. (b) A magnetic force microscopy 
image of a D= 1/im disk showing a magnetic vortex core at 
the center of the disk, (c) The polar Kerr signal as a function 
of pump-probe time delay measured for disks with diameters 
of 1.5 and 1 /im. (d) The vortex mode eigenfrequency as a 
function of the disk diameter. The analytical model (from 
Ref. Il3tl and simulation results are discussed in the text. 



the vortex deforms rapidly (i.e. on a time scale fast rel- 
ative to the core motion) in order to cancel out the edge 
charges 0. The experimental results are also in good 
agreement with a semi-analytical solution of the Landau- 
Lifshitz equation and micromagnetic simulations |l4l Il5| . 
The open circles in Fig. ^d) show the eigenfrequencies 
found in micromagnetic simulations of the full Landau- 
Lifshitz- Gilbert equation with a saturation magnetiza- 
tion M s = 770 emu/ cm 3 . damping parameter a = 0.008 
and a cell size of 5 nm[l6j. 

An additional feature of the data in Fig. ^c) is the 
higher frequency response at short timescales. As illus- 
trated in Fig.[2Ia), the torque M x H p has opposite signs 
in the two halves of a vortex. As a result, the spins in the 
upper half of the disk are rotated out of the plane, while 
those in the lower half of the disk are rotated into the 
plane. A spatially resolved polar Kerr image of a 1 fjjn 
diameter disk just after the pulse, demonstrating the ex- 
pected contrast, is shown in the right panel of Fig. |2Ja). 
The signal at short time-scales in Fig. ^c) is the subse- 
quent precession and decay of the non-equilibrium mag- 
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FIG. 2: (a) Schematic of the vortex magnetization (left), the 
torque exerted by the in-plane pulse (middle) , and polar Kerr 
image immediately after the pulse for a 1 ^m diameter disk 
(right), (b) The spectral power measured for four disk diame- 
ters with thicknesses of 20 nm. The curves are offset vertically 
for clarity. 



netization. The spectrum of these higher frequency spin 
wave modes can be obtained by Fourier transformation 
of the time-domain data. Spectra obtained for several 
different disk diameters are shown in Fig.[2Ib). 

Given the cylindrical symmetry of the vortex state, 
the spin wave modes can be indexed by the number 
of nodes n and m along the radial f and azimuthal 
4> directions: (n, m) = m„ jm (r;i)e'' m,1 " l "''"' m ^ where 
f(n,m) = Wn,m/27r is the eigenfrequency. This enumer- 
ation ignores the node that always exists at the vortex 
core, inside of which the in-plane components of the mag- 
netization vanish. The (0,0) mode is a purely radial mode 
and is excited most efficiently by an out-of-plane mag- 
netic field pulse as observed in Ref.|j. The (0,±1) modes 
are the two lowest order azimuthal modes and are excited 
by an in- plane pulsed fielddUGjJ. Because of the cylin- 
drical symmetry of the vortex ground state, the (0,±1) 
modes are expected to be degenerate if the core is fixed 
at the center of the disk. As expected, only a single peak 
is observed in the spectra of larger disks, as shown for 
the case of a 5 /Ltm diameter disk in Fig. EJb). 

Although a single (0,±1) mode is observed in the 5 /im 
diameter disk, the peaks in Fig.^b) begin to split as the 
disk diameter decreases. This doublet has also been ob- 
served in measurements on a 700 nm diameter disk by 
Zhu et aL0 The remainder of this paper is dedicated to 
understanding these doublets. Clearly the lifting of the 
degeneracy requires the breaking of cylindrical symme- 
try, and we argue that the reduced symmetry originates 
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FIG. 3: (color on-line) (a) The spectral power integrated over 
space for two different vortex-core polarities: blue (top) for 
the P = — 1 vortex state and red (bottom) for the P = +1 
state of a 1.5 /im diameter disk. The chiralities of the two 
vortex states are the same. The two curves are offset for clar- 
ity, (b) The spectral power and phase images of the (0,±1) 
modes at 3.2 GHz (top) and 3.8 GHz (bottom). The panels on 
the left show results for the P = +1 vortex state (power and 
phase, respectively) and the right panels the P = — 1 state. 
The arrows show the direction in which the phase changes 
from to 2ir. The response at 3.4 GHz (middle) is the super- 
position of the (0,-1) and (0,-1-1) modes. 



from the gyrotropic motion of the vortex core, an ar- 
gument supported by more rigorous calculations 0, Hsf 
and micromagnetic simulations To address this we 
will label the vortices by a polarity P = ±1, where P 
corresponds to the vortex core pointing out of {P = +1) 
or into (P = —1) the plane of the disk, and a chiral- 
ity C = ±1 indicating that the in-plane component of 
the magnetization rotates counterclockwise (C = +1) or 
clockwise (C = —1). 

As noted above, the low-frequency vortex mode is gy- 
rotropic and therefore has a particular sense of rotation. 
I gnor ing damping, the total force acting on the vortex 



F = G x 



dR dW(R) 



(1) 



dt <9R ' 
where R is the position of the core in the plane, G is 



the gyrovector with a direction (into or out of the plane) 
corresponding to the sign of the polarization P, and W is 
the total energy of the vortex. Since W does not depend 
on the chirality, the gyrotropic motion described by Eq.^ 
depends only the vortex polarity P. In our experiment, 
the core polarity can be controlled by relaxing the disk 
from saturation in the presence of a small (~ 80 Oe) out- 
of-plane field. The chirality cannot be controlled but is 
easily measured by imaging the displacement of the vor- 
tex in a small in-plane field. (The direction of displace- 
ment is perpendicular to the applied field and depends 
on the chirality.) In this way, we can select two vortices 
with identical chiralities but opposite polarities. Spectra 
for these two vortex states in a 1.5 /im disk are shown 
in Fig.[3Ja). Each spectrum shows three peaks at identi- 
cal frequencies of 0.2, 3.2, and 3.8 GHz. The upper two 
modes are the (0,±1) modes, for which a splitting can be 
resolved. 

Fig-EHb) shows the spectral power and phase images at 
three different frequencies for the two vortex polarities. 
These images are constructed from the Fourier transform 
of the time-domain datajj^]. As expected for azimuthal 
modes with one node, the phase wraps through 2-7T as the 
disk is traversed around the circumference. The spectral 
power images at the two peak frequencies (3.2 and 3.8 
GHz), shown in the top and bottom panels of Fig.|3fb), 
are essentially the same for the two polarities. The phase 
images, however, invert when the polarity of the vortex 
core is reversed. In other words, the direction in which 
the phase winds from — > 27r changes from counterclock- 
wise to clockwise. This corresponds to reversing the sign 
of the azimuthal index m, and so the (0,-1-1) and (0,-1) 
modes exchange identities when the polarity of the vortex 
core is reversed. 

It is also possible to examine the spectral power and 
phase in the frequency regime corresponding to a super- 
position of the two modes. This is shown in the middle 
panel of Fig. O^b). An equal superposition of the two 
(0,±1) modes would be a standing wave, and indeed a 
line of nodes can be seen in the spectral power image at 
3.4 GHz. Since the wave vectors of the (0,±1) modes are 
either parallel or anti-parallel to the magnetization, they 
have backward volume character (2^, and both eigenfre- 
quencies are therefore expected to lie below that of the 
radially symmetric (0,0) mode0,0|, although the (0,0) 
mode cannot be excited in the geometry of this experi- 
ment. 

As shown by the analysis of Fig. [3 the symmetry of the 
(0,±1) modes is determined by the polarity of the vortex 
core. The polarity determines the sense of rotation for 
the vortex gyrotropic motion, which breaks the azimuthal 
symmetry. Figure Ufa) shows the eigenfrequency of (0, 
±1) modes as a function of the dimensionless quantity 
(3 = L /R, where the thickness L is 20 nm and the radius 
R varies from 350 nm to 2.5 ^m. The solid and dashed 
curves show the frequencies of the two modes determined 
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FIG. 4: (a) The measured eigenfrequencies of the (0,±1) 
modes are shown as a function of /3 = L/R. The symbols 
are experimental data and the curves are the results of mi- 
cromagnetic simulations, (b) The splitting Af of the (0,±1) 
modes is shown as a function of the vortex mode frequency. 



theory all indicate that Af and fa are of the same order 
of magnitude, the discrepancies among the various ap- 
proaches are certainly significant. This reflects the dif- 
ficulty in treating exactly the coupling between the two 
different types of excitations. Since the coupling is mag- 
netostatic in origin, the experimental results may be sus- 
ceptible to non-ideal structure at the edges of the disk. 
This may also be addressed by carefully examining the 
role of the boundary conditions in the different theoreti- 
cal approaches. 

The broken degeneracy of the azimuthal spin-wave 
modes clearly demonstrates how low-frequency excita- 
tions associated with domain structure, of which a vortex 
is the simplest example, influence spin- wave dynamics. 
An intriguing question is how the coupling will evolve at 
larger wave vectors, corresponding to higher azimuthal 
eigenmodes. Since the magnetostatic interactions that 
provide the restoring force on the vortex core are long- 
range in nature, we expect that the coupling will be 
smaller for higher-order modes. At very large wave vec- 
tors, however, when the wavelength becomes comparable 
to the exchange length, spin-wave modes may interact 
strongly with the vortex core itself. 

This work was supported by the University of Min- 
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06029, the University of Minnesota Graduate School, and 
the Minnesota Supercomputing Institute. 



from micromagnctic simulations. As the radius increases, 
the eigenfrequencies of the (0, ±1) modes and the split- 
ting Af between the modes both decrease. For [3 < 0.01, 
A/ becomes smaller than the frequency resolution in our 
experiment and the (0,±1) modes are essentially degen- 
erate. Although both the absolute frequency and the size 
of the splitting are larger than the experimental values, 
the same trend is observed in the simulations. 

We have also examined the relationship between the 
magnitudes of the (0, ±1) splitting Af and the gyrotropic 
mode frequency fa- This is shown in Fig.^b) for both 
the experiment (closed circles) and simulations (open cir- 
cles). Over the range of our measurements, the splitting 
is approximately twice fa, although it is not possible to 
observe the gyrotropic mode in the largest disks. The 
simulations considered here as well as those of Zhu et 
ai.0| show a larger splitting relative to fa, although 
they also cannot be extended to the largest diameters, 
at which both Af and fa must approach zero. The 
relationship between the vortex core mode and the az- 
imuthal spin wave modes has also been considered re- 
cently by Ivanov and Zaspeldl. They calculate a split- 
ting Af which is approximately half fa, implying a curve 
that falls significantly below the experimental data in 
Fig. 0fb). Although experiment, simulations, and the 
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